Anisotropic microwave conductivity of cuprate superconductors in the presence of 

CuO chain induced impurities 
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The anisotropy in the microwave conductivity of the ortho-II YBa2Cu306.so is studied within the 
kinetic energy driven superconducting mechanism. The ortho-II YBa2Cu306.so is characterized by 
a periodic alternative of filled and empty b-axis CuO chains. By considering the CuO chain induced 
extended anisotropy impurity scattering, the main features of the anisotropy in the microwave 
conductivity of the ortho-II YBa2Cu3Oe.50 are reproduced based on the nodal approximation of 
the quasiparticle excitations and scattering processes, including the intensity and lineshape of the 
energy and temperature dependence of the a-axis and 6-axis microwave conductivities. Our results 
also confirm that the 6-axis CuO chain induced impurity is the main source of the anisotropy. 
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I. INTRODUCTION 



Impurity plays a crucial role in determining the 
behavior of many measurable properties of cuprate 
superconductors-^, such as the quasiparticle transport in 
the superconducting (SC) state. This follows from a fact 
that the physical properties of cuprate superconductors 
in the SC state are extreme sensitivity to the impurity 
effect than the conventional superconductors due to the 
finite angular-momentum charge carrier Cooper pairing 
with the d-wave symmetry-i. The single common fea- 
ture is the presence of the Cu02 plane- and it seems 
evident that the unusual behaviors of cuprate supercon- 
ductors are dominated by this Cu02 planed. Experimen- 
tally, over twenty years measurements of electrodynamic 
properties at microwave energies have provided rather 
detailed information on the quasiparticle transport of 
cuprate superconductor o 4 ' 5 ' 6 ' 7 ' 8 ' i, where the nodal quasi- 
particle spectrum contains most of the features expected 
for weak-limit impurity scattering. The early microwave 
conductivity measurements 4 **^ 7 - showed the microwave 
conductivity spectrum is Cu02 plane isotropic. However, 
the recent development of a sufficient number of fixed- 
energy cavity perturbation system to map out a coarse 
microwave conductivity spectrum^ allowed to resolve 
additional features in the microwave conductivity spec- 
trum. Among these new achievements is the observation 
of an anisotropy in the microwave conductivity spectrum 
of the highly ordered ortho-II YP^CusOVso, where the 
well pronounced difference in intensity and lineshape in 
the a-axis and 6-axis directions is remarkable. 

Although the anisotropy in the microwave conduc- 
tivity spectrum of the ortho-II YBa2Cu306.5o is well- 
established experimentally^, its full understanding is 
still a challenging issue. Theoretically, based on a phe- 
nomenological Bardeen-Cooper-Schrieffer (BCS) formal- 
ism with the d-wave SC gap function, it has been ar- 
gued in the mean-field level that the interband tran- 
sitions produce a strongly anisotropic feature 10 . How- 
ever, it has been shown that the Cu02 bilayer may lead 
to two nearly identical contributions to both a-axis and 
6-axis microwave conductivities*-, since the splitting be- 
tween bonding and antibonding combinations of the pla- 
nar wavefunctions is expected to be weak- 11 . In our 



earlier worki*- based on the kinetic energy driven SC 
mechanis m 13 i 14 , the effect of the extended impurity scat- 
tering potential on the quasiparticle transport in the non- 
ortho-II phase of cuprate superconductors has been dis- 
cussed within the nodal approximation of the quasiparti- 
cle excitations and scattering processes, and the obtained 
energy and temperature dependence of the microwave 
conductivity are consistent with the experimental data in 
the non-ortho-II phase of cuprate superconductors 4 ^^ 7 -. 
In the d-wave SC state of cuprate superconductors, the 
characteristic feature is the existence of four nodal points 
[±7r/2, ±7r/2] (in units of inverse lattice constant) in the 
Brillouin zone2iI£, where the SC gap function vanishes, 
therefore the Cu02 plane currents are mainly carried 
by nodal quasiparticles, and the quasiparticle transport 
properties of cuprate superconductors in the SC state are 
largely governed by the quasiparticle excitations around 
the nodes. Since the gap nodes lie close to the diago- 
nals of the Brillouin zone, these quasiparticle excitations 
carry both a-axis and 6-axis currents, then the microwave 
conductivity should be Cu02 plane isotropic. However, 
the ortho-II YBa 2 Cu 3 06.5o is a very special cuprate SC 
material, and is characterized by a periodic alternative of 
filled and empty 6-axis CuO chains*^. In particular, the 
recent experimental measurements^ unambiguously es- 
tablish that the 6-axis CuO chain induced impurity is the 
dominant source of the anisotropy in the microwave con- 
ductivity spectrum of the ortho-II YBa2Cu306.so- This 
experimental result also implies that the 6-axis CuO 
chain in the ortho-II YBa2Cu30g.5o induced impurities 
lead to an anisotropy of the extended impurity scat- 
tering potential in the Cu02 planes. In this paper we 
show explicitly if the effect of the anisotropy of the ex- 
tended impurity scattering potential induced by the 6- 
axis CuO chain is considered within the kinetic energy 
driven SC mechanism, one can reproduce some main 
features of the anisotropy in the microwave conductiv- 
ity spectrum observed*^ experimentally on the ortho-II 
YBa 2 Cu 3 6 .5o. 

The rest of this paper is organized as follows. In 
Sec. II, we present the basic formalism, where the BSC- 
like Green's function under the kinetic energy driven SC 
mechanism is dressed via the extended anisotropy im- 
purity scattering. Within this framework, we calculate 
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explicitly the a-axis and 6-axis microwave conductivi- 
ties based on the nodal approximation of the quasipar- 
ticle excitations and scattering processes. The energy 
and temperature dependence of the a-axis and o-axis mi- 
crowave conductivities of the ortho-II YBa2Cu 3 Og.5o are 
presented in Sec. III. Sec. IV is devoted to a summary. 



II. THEORETICAL FRAMEWORK 

The basic element of cuprate superconductors is two- 
dimensional CuC>2 planes^ as mentioned above. It has 
been shown that the essential physics of the doped Cu02 
plane is properly accounted by the t-J model on a square 
lattice 2 -^, 



H = -* E ct^w + a* E °l c ™ + 3 E s » • s 



(1) 



where i) — ±i, ±y, C| (G^) is the electron creation (an- 
nihilation) operator, = (Sf , Sf , Sf) is spin operator, 
and /i is the chemical potential. This t-J model is sub- 
ject to an important local constraint ^ C ia Ci a < 1 to 
avoid the double occupancy 3 -, which can be treated prop- 
erly in analytical calculations within the charge-spin sep- 
aration (CSS) fermion-spin theor y 14 ' 16 , where the con- 
strained electron operators are decoupled as = S~ 
and Cii = h\,Sf, with the spinful fermion operator 

h ia = e^^hi describes the charge degree of freedom to- 
gether with some effects of spin configuration rearrange- 
ments due to the presence of the doped charge carrier 
itself, while the spin operator S% describes the spin de- 
gree of freedom, then the electron local constraint for the 
single occupancy is satisfied in analytical calculations. In 
this CSS fermion-spin representation, the t-J model (1) 
can be expressed as, 

if) 

- M h\ a h ia + Jeff ^ Si ■ Sj+fj, (2) 



with J e g = (1 — 5) 2 J, and 5 = {hl a h icr ) — {h\hi) is the 
charge carrier doping concentration. For a understand- 
ing of the SC state properties of cuprate superconduc- 
tors, the kinetic energy driven SC mechanism has been 
develope d 13 ' 14 , where the interaction between charge car- 
riers and spins from the kinetic energy term in the t- 
J model (2) induces the charge carrier pairing state 
with the d-wave symmetry by exchanging spin excita- 
tions, then the electron Cooper pairs originating from 
the charge carrier pairing state are due to the charge- 
spin recombination, and their condensation reveals the 
SC ground-state. In particular, it has been shown that 
this SC state is a conventional BCS like with the d-wave 
symmetry^, so that the basic BCS formalism with the 
d-wave SC gap function is still valid in quantitatively re- 
producing all main low energy features of the SC coher- 
ence of quasiparticles, although the pairing mechanism is 
driven by the kinetic energy by exchanging spin excita- 
tions. Following our previous discussions—, the electron 
Green's function in the SC state can be obtained in the 
Nambu representation as, 



G(k,w) 



ujtq + A z Qt)n + g k r 3 
to 2 -El 



(3) 



where To is the unit matrix, t\ and t 3 are Pauli matrices, 
other notations are defined as same as in RefJ^, and have 
been determined by the self-consistent calculatio n 13 ' 14 . 

In the presence of impurities, the unperturbed elec- 
tron Green's function (3) is dressed via the impurity 
scatterin g 18 ' 19 ' 20 , 

G/(k,w) = 22Gi a T a (k,uj) 

= [Gfow)- 1 -^.*)]- 1 , (4) 

with the self-energy E(k, u>) — E a (k, w)r a . It has 
been shown that all but the scalar component of the 
self-energy function can be neglected or absorbed into 
Az(k ) 18 ' 19 ' 20 , In this case, the dressed electron Green's 
function (4) can be explicitly rewritten as, 



[w - E (k, u)]t + A z (k)n + [si. + S 3 (k,cj)]T 3 



Eo(k,w)]a 
I 



A|(k) 



(5) 



where the self-energies Eo(k, cu) and X 3 (k, lu) are treated 
within the framework of the T-matrix approximation as, 



E(k,w) = /OiT kk (w), 



(6) 



where pi is the impurity concentration, and Tkk(w) is the 
diagonal element of the T-matrix, 

Tkk'M = Vkk'73 + J2 ^kk"T 3 G/(k", w)£w(w), (7) 

k" 

with Vkk' is the impurity scattering potential and 
Ikk'(w) = T^ k ,(cj)r + T^ k ,(w)r 3 . In our earlier 



work without considering the &-axis CuO chain induced 
impurities^ 2 -, we have discussed the effect of the extended 
isotropy impurity scattering potential on the quasiparti- 
cle transport in the non-ortho-II phase of cuprate super- 
conductors within the nodal approximation of quasiparti- 
cle excitations and scattering processes, where there is no 
gap to the quasiparticle excitations at the four nodes, and 
then the quasiparticles are generated only around these 
four nodes. In this CUO2 plane isotropic case, a gen- 
eral scattering potential Vkk' need only be evaluated in 
three possible case o 18 ' 19 : the intranode impurity scatter- 
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FIG. 1: The anisotropic impurity scattering within the d- 
wave case. Vi, V2, V3, and Vi are potentials for intranode, 
adjacent-node along the a-axis direction, opposite-node, and 
adjacent-node impurity scattering along the b-axis direction 
scattering. 



quasiparticles moving in the CuC>2 planes^. After incor- 
porating this 6-axis CuO chain induced anisotropy impu- 
rity scattering into the extended isotropy impurity scat- 
tering potential, the total scattering potential Vkk' need 
be evaluated in four possible cases as shown in Fig. 1: 
the intranode impurity scattering Vkk' = V\ (k and k' at 
the same node) and the opposite-node impurity scatter- 
ing Vkk' = V3 (k and k' at the opposite nodes), these two 
cases are the same as in the previous discussions in the 
non-ortho-II phase of cuprate superconductors 12 . How- 
ever, the adjacent-node impurity scattering along the a- 
axis direction Vkk' = V2 (k and k' at the a-axis adjacent 
nodes) is different from that along the b- axis direction 
Vkk' = V4 (k, and k' at the 6-axis adjacent nodes) in 
the present discussions of the ortho-II YBa2Cu306.5o- In 
this case, the impurity scattering potential Vkk' in the 
T-matrix is effectively reduced as, 



ing Vkk' = Vi (k and k' at the same node), the adjacent- 
node impurity scattering Vkk' = V2 (k and k' at the ad- 
jacent nodes), and the opposite- node impurity scattering 
Vkk' = V3 (k and k' at the opposite nodes). However, 
this CuC>2 plane isotropic case is broken in the presence 
of the 6-axis CuO chain induced impurities in the ortho- 
II YBa2Cu306.5o, since the CuO chain induced impurity 
scattering potential can be described by an anisotropic 
potential in the Cu02 planes^. This simply assumes 
that the screened Coulomb potential created by the CuO 
chain induced impurities produces a footprint sensed by 



Vkk' -> V = 



Vi V 2 V 3 V 4 

v 2 vi v 4 v 3 

V 3 V 4 Vi V2 

V A V 3 V 2 Vi 



(8) 



We emphasize that the anisotropy of the impurity scat- 
tering potential (then Cu02 plane microwave conductiv- 
ity) has been reflected by this important difference be- 
tween the adjacent-node impurity scattering strengths V 2 
and V4. Now we follow our previous discussions^, and 
obtain explicitly the a-axis and 6-axis microwave conduc- 
tivities of the ortho-II YE^CusOe+y as, 



cr a (w,T) = 4e 2 i; 2 . / P J -[ReJ a (uj' - iS, J + oj + iS) - ReJJa/ + iS, J + w + iS) ], (9a) 

J-00 2?r w 

a b (cj,T) = Ae 2 v 2 f [°° — ~ "f^' + w) [ReJ b (a/ -iS,u' + u} + iS) - ReJ 6 (a/ + i5,ui' + u} + iS)], (9b) 

J-00 2tt u' 



where tif{lo) is the fermion distribution function, Vf — kernel function J a (uj' ,uj) and J b (u>' ,u>) are given by, 
spit is the electron velocity at the nodal points, and the 



Ja(uJ ' ,UJ) = [i - (l«4 0) + Lg/f )][i - fef + lyb)] - Vk ] + l«A 3) wA 0) + v¥] ' (10a) 

j(0) + L b rj(0) j(3) + j(3) j(0h 

Jb{LJ ' ,uj) = [i - we + - (Vf + hh] ~Wi^ + L\if] [L \^ + l»4 0) ] ' (10b) 



r 

where the functions LI, L%, L\, and if are expressed as, 



- T? 3 (u>')T? 3 {u' + u) - TU^)TUu' + cj) + TlWU"' +")+ T? 4 ("')T? 4 ("' + (Ha) 
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u>) = pi [T X1 (u)')Tf x (w' + + I* (w')!*^' + cj) - T? 2 (lj')T? 2 (lj> + u) - Tf 2 (u/)T° 2 (a/ + u) 

- T? 3 (u>')T? 3 (u' + W ) - T^u/)^' + W ) + I* (a;')^' + w) + T^T^cj' + u,)], (lib) 

L\ (w', w) = pi [T^ (ci/JT?! (w' + w) + 7 1 ! 3 ! (a/)??! (w' + tw) + T? 2 ( W ')^12 + w) + I? a (w')2u (w' + w) 

- T^u/)^' + w) - TUu')TUu' + u) - T^{J)Tl{J + u) - TUu')TUu' + (He) 
^2 w) = ft Pft (w')Tj\ {of + oj) + Tl (co')T^ (lo> + u ) + T° 2 (u/)T& (w' + w) + T x 3 2 ( w ')3? 2 (w' + w) 

- ^3(^)^3^' + w) - T^u/)^' + u) - T° 4 (u;')TUu;' TU^TU^' + ")], (Hd) 



and the functions Jg (w, a/ )and4 u >,w') are evaluated in terms of the dressed Green's function (5) as, 



C(uj,uj')t + 4 3 \tu,u J ')T 3 = -J2Gi0t^)G I (KLo + Lu>), (12a) 

k 

lf\u,u')r + 4 3) ( w , w ')T3 = ^^G/(k^)f3G / (k^ + c i ;'). (12b) 



It is clearly that if the effect of the 6-axis CuO chain 
induced impurities is neglected, i.e., V 2 = V4, this leads 
to LI(lo',uS) = L\{u)',uj) and L%{w , ,w) = L b 2 (uj',u}), and 
then the d-axis and 6-axis microwave conductivities in 
Eq. (9) are reduced to the isotropic onei2 a a (u)',ui) = 
<jb(u)',ui) — <t(lu' ,u), 

III. ENERGY AND TEMPERATURE 
DEPENDENCE OF THE fi-AXIS AND 6-AXIS 
MICROWAVE CONDUCTIVITIES FOR THE 
ORTHO-II YBA2CU3O6.50 

In cuprate superconductors, although the values of J 
and t is believed to vary somewhat from compound to 
compound 2 , however, as a qualitative discussion, the 
commonly used parameters in this paper are chosen 
as t/J = 2.5, with an reasonably estimative value of 
J ~ 1000K. We are now ready to discuss the energy and 
temperature dependence of the d-axis and 6-axis quasi- 
particle transport of the ortho-II YBa 2 Cu306.5o with the 
extended anisotropy impurity scattering. We have per- 
formed a calculation for the energy dependence of the 
d-axis and 6-axis microwave conductivities cr a (w,T) and 
<7b{oJ,T) in Eq. (9) at low temperatures, and the results 
of a a (u>,T) (top panel) and <7&(w,T) (bottom panel) as 
a function of energy with temperature T = 0.002 J = 2K 
(solid line), T = 0.004 J = 4K (dash-dotted line) and 
T = 0.008 J = 8K (dashed line) under the slightly strong 
impurity scattering potential with V\ = 100 J, V 2 = 90 J, 
V3 = 30J, and V4 = 30J at the impurity concentration 
Pi = 0.00002 for the doping concentration 5 = 0.15 are 
plotted in Fig. 2 in comparison with the correspond- 
ing experimental results^ of the ortho-II YBa 2 Cu306.5o 
(inset). It is clearly that the anisotropy of the en- 
ergy evolution of the low temperature microwave con- 
ductivity of the ortho-II YBa 2 Cu306.so is qualitatively 
reproduced®^., where quasiparticle spectral weights dif- 
fer by a factor of two between the d-axis and 6-axis di- 



rections. Moreover, although the cusplike lineshapes are 
observed along both d-axis and 6-axis directions as in 
the previous case in the non-ortho-II phase of cuprate 
superconductors^ 2 -, the width of the d-axis spectrum is 
significantly broadened. This broadening is only at- 
tributed to increased quasiparticle scattering arising from 
the CuO chain induced impurities. In comparison with 
our previous results for the non-ortho-II phase of cuprate 
superconductors 1 ^, the present results of the anisotropy 
therefore confirm that the CuO chain induced impurity is 
the dominant source of the anisotropy in the microwave 
conductivity spectrum of the ortho-II YBa 2 Cu306.5o 9 - 

In the above discussions of the low temperature and 
low energy case, the inelastic quasiparticle-quasiparticle 
scattering process has been dropped, since it is sup- 
pressed at low temperatures and low energies due to the 
large SC gap parameter in the quasiparticle excitation 
spectrum. For a better understanding of the anisotropy 
in the microwave conductivity spectrum, we now turn to 
discuss the temperature dependence of the quasiparticle 
transport of the ortho-II YBa 2 Cu306.5o, where T may 
approach to T c from low temperature side, and therefore 
the inelastic quasiparticle-quasiparticle scattering pro- 
cess should be considered 1 ^. The contribution from this 
inelastic quasiparticle-quasiparticle scattering process is 
increased rapidly when T approaches to T c from low tem- 
perature side, since there is a small SC gap parameter 
near T c . In particular, it has been pointed outj 2 ^ that the 
contribution from the quasiparticle-quasiparticle scatter- 
ing process to the transport lifetime is exponentially sup- 
pressed at low temperatures, and therefore the effect 
of this inelastic quasiparticle-quasiparticle scattering can 
be considered by adding the inverse transport lifetime 20 
TjZzAT) to the imaginary part of the self-energy function 
Hq(uj) in Eq. (5) as in our previous discussions^ 2 -, then 
the total self-energy function £ ot (u;) can be expressed 
asi 2 -^ E tot (w) = Sq(lj) _ i[ 2Tincl (T)}-\ with T inol (T) 
has been chosen as prme^T)]- 1 = 2 x 10 4 (T - 0.002) 4 J. 
Using this total self-energy function £ ot (u>) to replace 
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FIG. 2: The microwave conductivities of the a-axis (top 
panel) and 6-axis (bottom panel) as a function of energy with 
T = 0.002 J = 2K (solid line), T = 0.004 J = 4K (dash-dotted 
line) and T = 0.008J = 8K (dashed line) at p; = 0.00002 for 
t/J = 2.5 with Vi = 100J, V 2 = 90 J, V 3 = 30J, and V 4 = 30 J 
in S = 0.15. Inset: the corresponding experimental result of 
the ortho-II YBa2Cu3C>6.50 taken from Ref. H|. 



So 0*0 i n Eq. (5), we have performed a calculation for 
the temperature dependence of the a-axis and o-axis mi- 
crowave conductivities a a (uj,T) and <7b(w,T), and the 
results of a a {to,T) and <7b(w,T) as a function of tem- 
perature T with energy to = 0.0000547J w 1.14GHz 
(solid line), w = 0.0001094 J w 2.28GHz (dash-dotted 
line), w = 0.0006564 J w 13.4GHz (dashed line), and 
w = 0.001094J w 22.8GHz (dotted line) under the 
slightly strong impurity scattering potential with Vi = 
100J, V 2 = 90J, = 30J and V 4 = 30J at the im- 
purity concentration pi = 0.00002 for the doping con- 
centration S = 0.15 are plotted in Fig. 3. For com- 
parison, the corresponding experimental results- of the 
ortho-II YBa2Cu306.5o are also plotted in Fig. 3 (in- 
set). Obviously, the overall temperature dependence is 
similar to that of the temperature dependence of the mi- 
crowave conductivity in the non-ortho-II phase of cuprate 
superconductors^, where both a-axis and o-axis temper- 
ature dependent microwave conductivities a a (uj,T) and 
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FIG. 3: The microwave conductivities of the a-axis (top 
panel) and 6-axis (bottom panel) as a function of temperature 
with energy u » 1.14GHz (solid line), uj ~ 2.28GHz (dash- 
dotted line), to w 13.4GHz (dashed line), and w « 22.8GHz 
(dotted line) at p l = 0.00002 for t/J = 2.5 with Vi = 100J, 
V 2 = 90J, V 3 = 30J, and Va = 30J in 8 = 0.15. Inset: the cor- 
responding experimental result of the ortho-II YBa2Cu306.so 
taken from Ref. 8. 



ab(u>,T) increases rapidly with increasing temperatures 
to a broad peak, and then falls roughly linearly. In par- 
ticular, this broad peak shifts to higher temperatures as 
the energy increased. However, the spectral width is evi- 
dently anisotropic, where the spectral width of the d-axis 
microwave conductivity cr a (a->,T) is larger than the cor- 
responding value of the o-axis one. Moreover, the peak 
position of the a-axis microwave conductivity a a (u>,T) 
is located at higher temperature than the corresponding 
value of the 6-axis one, in qualitative agreement with the 
experimental data of the ortho-II YBa2Cu306. sere- 
in our present theory, the anisotropy of the microwave 
conductivity reflects directly from the electron vertex cor- 
rection due to the extended anisotropy impurity scatter- 
ing potential. However, it has been shown that for all 
scattering strength the thermal vertex correction is neg- 
ligible compared to the electron onei&, therefore we ex- 
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pects the heat transport of the ortho-II YBa2Cu30g.5o to 
be C11O2 plane isotropic. This is a unique feature of our 
theory that is different from the other theory^ in which 
the fermi surface is assumed to be altered and thus the 
heat transport may be also anisotropic. This should be 
verified by further experiments. 

IV. SUMMARY 

Within the framework of the kinetic energy driven SC 
mechanism, we have studied the anisotropy in the mi- 
crowave conductivity spectrum recently observed^ 9 - in 
the ortho-II YBa 2 Cu 3 06.5o- The ortho-II YBa 2 Cu 3 6 .5o 
is characterized by a periodic alternative of filled and 
empty &-axis CuO chains. The extended anisotropy im- 
purity scattering potential results when the extended 
isotropy impurity scattering potential in the Cu0 2 planes 
incorporates the 6-axis CuO chain induced impurity scat- 
tering potential. Based on the nodal approximation of 



the quasiparticle excitations and scattering processes, we 
have calculated the a-axis and 6-axis microwave conduc- 
tivities with the extended anisotropy impurity scattering, 
and reproduced the main features^ of the anisotropy 
in the microwave conductivity spectrum of the ortho-II 
YBa 2 Cu306.5o, including the intensity and lineshape of 
the energy and temperature dependence of the a-axis and 
6-axis microwave conductivities. Our results also confirm 
that the 6-axis CuO chain induced impurity is the main 
source of the anisotropy 9 . 



Acknowledgments 

This work was supported by the National Natural Sci- 
ence Foundation of China under Grant No. 10774015, 
and the funds from the Ministry of Science and Tech- 
nology of China under Grant Nos. 2006CB601002 and 
2006CB921300. 



To whom correspondence should be addressed. 

1 See, e.g., the review, H. Alkral, J. Bobroff, M. Gabay, and 
P. J. Hirschfeld, Rev. Mod. Phys. 81, 45 (2009), and ref- 
erences therein. 

2 See, e.g., the review, A. Damascelli, Z. Hussain, and Z. 
X. Shen, Rev. Mod. Phys. 75, 473 (2003), and references 
therein. 

3 P. W. Anderson, in Frontiers and Borderlines in Many 
Particle Physics, edited by R. A. Broglia and J. R. Schrief- 
fer (North-Holland, Amsterdam, 1987), p. 1; Science 235, 
1196 (1987). 

4 D. A. Bonn, R. Liang, T. M. Riseman, D. J. Baar, D. C. 
Morgan, K. Zhang, P. Dosanjh, T. L. Duty, A. MacFarlane, 
G. D. Morris, J. H. Brewer, W. N. Hardy, C. Kallin, and 
A. J. Berlinsky, Phys. Rev. B 47, 11314 (1993). 

5 Shih-Fu Lee, D. C. Morgan, R. J. Ormeno, D. M. Broun, 
R. A. Doyle, J. R. Waldram, and K. Kadowaki, Phys. Rev. 
Lett. 77, 735 (1996). 

6 A. Hosseini, R. Harris, Saeid Kamal, P. Dosanjh, J. Pre- 
ston, Ruixing Liang, W. N. Hardy, and D. A. Bonn, Phys. 
Rev. B 60, 1349 (1999). 

7 P. J. Turner, R. Harris, Saeid Kamal, M. E. Hayden, D. 
M. Broun, D. C. Morgan, A. Hosseini, P. Dosanjh, G. K. 
Mullins, J. S. Preston, Ruixing Liang, D. A. Bonn, and W. 
N. Hardy, Phys. Rev. Lett. 90, 237005 (2003). 

8 R. Harris, P. J. Turner, Saeid Kamal, A. R. Hosseini, P. 
Dosanjh, G. K. Mullins, J. S. Bobowski, C. P. Bidinosti, D. 
M. Broun, Ruixing Liang, W. N. Hardy, and D. A. Bonn, 
Phys. Rev. B 74, 104508 (2006). 

9 J. S. Bobowski, P. J. Turner, R. Harris, Ruixing Liang, 



D. A. Bonn, and W. N. Hardy, Physica C 460-462, 914 
(2007) ; |arXiv:cond-mat/0612344| 

E. Bascones, T. M. Rice, A. O. Shorikov, A. V. Lukoyanov, 
and V. I. Anisimov, Phys. Rev. B 71, 012505 (2005). 

1 W. A. Atkison, Phys. Rev. B 59, 3377 (1999); Yu Lan, 
Jihong Qin, and Shiping Feng, Phys. Rev. B 76, 014533 
(2007). 

2 Zhi Wang, Huaiming Guo, and Shiping Feng, Physica C 
468, 1078 (2008). 

3 Shiping Feng, Phys. Rev. B 68, 184501 (2003); Shiping 
Feng, Tianxing Ma, and Huaiming Guo, Physica C 436, 
14 (2006). 

4 See, e.g., the review, Shiping Feng, Huaiming Guo, Yu Lan, 
and Li Cheng, Int. J. Mod. Phys. B 22, 3757 (2008). 
See, e.g., the review, C. C. Tsuei and J. R. Kirtley, Rev. 
Mod. Phys. 72, 969 (2000). 

6 Shiping Feng, Jihong Qin, and Tianxing Ma, J. Phys. Con- 
dens. Matter 16, 343 (2004). 

7 Huaiming Guo and Shiping Feng, Phys. Lett. A 361, 382 
(2007). 

8 A. C. Durst and P. A. Lee, Phys. Rev. B 62, 1270 (2000). 

9 Tamara S. Nunner and P. J. Hirschfeld, Phys. Rev. B 72, 
014514 (2005). 

D. Duffy, P. J. Hirschfeld, and D. J. Scalapino, Phys. Rev. 
B 64, 224522 (2001). 

1 S. Graser, P. J. Hirschfeld, and L. Y. Zhu, Phys. Rev. B 
76, 054516 (2007). 

2 M. B. Walker and M. F. Smith, Phys. Rev. B 61, 11285 
(2000). 



